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Abstract 

We discuss the gravitational dual of a holographic superconductor consisting of a U(l) gauge field, 
a complex scalar field coupled to a charged AdS black hole and a higher-derivative coupling between 
the U(l) gauge field and the scalar with coupling constant n. In the presence of a magnetic field, 
the system possesses localized spatially dependent droplet solutions which, in the low temperature 
limit, have smaller critical temperature for r\ > than the droplet solutions without the interaction 
term (n = 0). In the weak magnetic field limit, the opposite behavior is observed: the critical 
temperature increases as we increase r\. We also calculate the energy gap in the probe limit and 
find that it is larger for n < than the energy gap in the conventional case (n = 0). 

PACS numbers: 11.25.Tq, 04.70.Bw, 74.20.-z 
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I. INTRODUCTION 



AO 



The plicnomonological holographic superconductor [1|, |2| has many similarities with the effective Landau- 
Ginzburg theory. The gravity sector of the holographic superconductor is an Abelian-Higgs model with a 
dynamical black hole metric and the scalar field condenses below a certain critical temperature. The Landau- 
Ginzburg theory is the effective field theory description, in one less dimension, of a superconductor near the 
superconducting phase transition. The dynamical degree of freedom is the complex order parameter, ^ = (0), 
which is coupled to a background electromagnetic field. This is described by the standard Ginzburg-Landau 
functional 

F = a\^ + ^\\ (1) 

and the coefficient a vanishes at the transition temperature T = T c . At T < T c , the coefficient a is negative 
and the minimum of F occurs for a uniform superconducting state at |^| 2 = — ^. 

However, the Landau-Ginzburg theory is not a microscopic theory and for the description of the dynamics 
of a Cooper pair the BCS theory is required. By the same token, one hopes that the phenomenological 
description of holographic superconductors will find its microscopic realization in string theory 3|, m- I n 
spite of this limitation, if one insists on working with the Landau-Ginzburg theory and wishes to describe 
superconductivity away from the critical temperature and in the presence of a magnetic field.the functional 
([T]) has to be supplemented with higher order terms involving spatial derivatives of the form 5| 

Fq = a(B, T)\V\ 2 + 7 (g, T)|^| 2 + ^ T) |^ 2 fr| 2 + ^ T) |^| 4 . (2) 

Note that the various coefficients also depend on the magnetic field. The critical temperature resulting from 
the generalized functional (|2|) depends on the value of the coefficient 7 and for 7 < interesting new physics 
emerges 6 1 . 

Despite the similarities between the holographic superconductors and the Landau-Ginzburg theory, there 
are important differences discussed in [2(. However, to keep the similarities of the holographic superconductor 
with the generalized Landau-Ginzburg theory, we can add a potential term which depends on the scalar field 
'J, as well as its spatial derivatives along the boundary, 

y(*,^ 2 *,V 2 *,...) . (3) 



While spatial derivatives are used in @, they are problematic in our case, because they signal the breakdown 
of Lorentz invariance. To obtain an expression similar to @ in a covariant way, we note that the vector 
F^ v d u ^ consists of the derivatives of ^> along the x and y directions, if the electromagnetic field strength F M „ 
describes a magnetic field in the radial direction (perpendicular to the boundary) . This leads us to consider 



(4) 



where F = dA is the field strength of a U(l) gauge field A^, ^ is a charged complex scalar field of charge q 
and mass m and = — iqA^. 

In this work we study the effects of the presence of the higher-derivative terms on the holographic super- 
conductor. To avoid technical difficulties, we set all couplings, except r\ to zero (rf = ■ ■ ■ = 0). In section 
HTl we present the details of the model. In section IIII1 we consider droplet solutions and calculate the critical 
temperature. In section IIV1 we go below the critical temperature and calculate the energy gap in the probe 
limit. Finally in section |V1 we discuss our conclusions. 



II. THE MODEL 



We consider the action 



S = d 4 c 



R + 6/L 2 1 



-F^F^ - \D^\ 2 - VpS!) 
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where the potential term is given by ^ with rj' = ■ ■ ■ = 0. The field equations are: 
• the Einstein equations 



R 



1 3 



(5) 



(6) 



• the Maxwell equations 



V^F^ + ^=<9 M [y/^g{D K ^f)(D x ^y (g KV F» x - g K »F uX + g^F^ - g^ x F UK )] 
= iq - + iqr]g iip F pv [F^*(D K ^) - F^ A *(L> A *)*] 



(7) 



• and the scalar field equation 
1 
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d ll [y/=$g> iv {d v V - iqA v *)] + iqgT ' A v (d^ - iqA^) + m 2 4- 



V 



■A [y/=g~g KX F KV F x » {d„* - iqA v ^)] + iqr,g KX F KV F x v A„ (d v V - iqA,^) = . (8) 



In this work, we will set L = 1, 8irG = 1, g = 1. 

We note that the presence of the coupling constant r\ adds new terms in the field equations which makes 
the system of the differential equations highly non-trivial. However, in spite of the complexity of the system, 
we were able to perform numerical calculations and also calculate the critical temperature analytically in the 
weak field limit. 

III. THE CRITICAL TEMPERATURE 



The original motivation of introducing a potential term of the form ^ in the Landau-Ginzgurg theory was 
to study the behavior of a superconductor in the presence of a strong magnetic field and the generation of 
anisotropic states in the low temperature limit. To investigate further similarities of a holographic supercon- 
ductor to the generalized Landau-Ginzburg theory, in this section we switch on a magnetic field and study 
the effect of spatially dependent solutions on the critical temperature. 



Spatially dependent solutions in the presence of a magnetic field were studied in [9h11| . Solving the Maxwell 



scalar equations in a dyonic black hole background localized droplet and vortex solutions were found and their 



(B,T) phase diagram was studied. In 



with the addition of an external magnetic field, it was argued 



that the holographic superconductor is of type II as the magnetic field is lowered. An analytic study of a 



holographic superconductor was carried out in 



12j . where the upper critical magnetic field was calculated 



analytically. Vortex lattice solutions in the presence of a magnetic field were presented in 13 1. 

To study the effect of the interaction term (with r\ ^ 0), we will consider the backreaction of the Maxwell 
field to gravity. The Einstein-Maxwell field equations ((6]) and (0 at the critical temperature are solved by a 



dyonic black hole (i.e., one with both an electric and a magnetic charge). The metric reads 



14] 



ds 2 



dr 2 



+ r 2 [-h{r)dt 2 + dx 2 + dy 2 



h(r] 



(X 2 + B 2 )r% (I 
4r 3 



The gauge field is given by 



Aft = A t {r)dt + A y (x)dy = Ar+ ( 1 - ^ ) dt + Br\xdy , 



(9) 



(10) 

where according to the AdS / CFT dictionary, the charge density and the chemical potential on the boundary 
are given by p = Ar^_ and /x = Ar + , respectively. Br 2 ^ is the magnitude of the external magnetic field (in the 
radial direction, r, perpendicular to the boundary xy-plane). At the horizon, r — r+, we have At = 0. The 
Hawking temperature of the black hole is 



3r + / A 2 +B 2 
4tt I 12 



(11) 



In the above background with a magnetic field B, we will investigate the spatially dependent solutions of the 
field equations by considering the ansatz ^ = ^(r, x). Near the AdS boundary, the asymptotic behavior is 
^ = -^L- + with A+ = | ± i-v/9 + 4m 2 . Both and "f? correspond to normalizablc modes, therefore we 

r — T ' J, Z 

obtain two different systems (labeled by A±) in which one mode is the source and the other gives the vacuum 
expectation value of the dual operator of dimension Aj ('J'i ~ (Cj))- 

Right below the critical temperature, we may regard the scalar field as a perturbation in the gravitational 
background of © and (|10l) . Then the scalar field equation in the coordinate z 



1 + 

h 













-!)■ 







,4 /i/2 iJ2 



becomes 



(1 V z*A?/t2)Al gV ^ v 



ft 



(12) 



where the prime denotes differentiation with respect to z. Without loss of generality, we can take to be real 
and separate 'J into functions of a single variable, 



Combining equations (fT2]) and (p~3|) . we obtain 
ft 



(13) 



ft' 2 
"ft ~ "z 



ft*; 



(i-')^; 2 K)^ 



?*ift 



^(l+^V) 

*2 L J 



For consistency of this equation, we need 



#2 - £ 2 z 2 *2 = -fc 2 ^ 



(14) 



(15) 



where k 2 is an arbitrary parameter. This equation is of the form of Schrodinger's equation for a simple 
harmonic oscillator. The eigenvalues are 



fc 2 = (2n+ l)B , n = 0,1,2,. 



(16) 



and the corresponding eigenfunctions are given in terms of Hermite polynomials, H n (y2Bx). Note that the 
x-dependent part of the scalar field obeys the same equation as the one found in 9|. As explained in {J, we 
ought to choose the ground state (n = 0) to find the critical temperature. 

Substituting the x-dependent equation (JT5J) into (|T4)) . we obtain the field equation for ty\{z) which reads 

'ft' 2 



[1 - V X 2 z 4 ] *i' 
2 

j2ft 



2 „4 



r/X z 



ft' 2 

ft z y ' \ h z 
X 2 (l-i 1 X 2 z 4 )(l-z) 2 _ 5(1 + ?? ff 2 z 4 ) 
ft 2 ft 



*i = , 



(17) 



FIG. 1: The ^-dependence of the scalar field for various values of the magnetic field B and coupling constant rj. 

where we have chosen m 2 = — 2 for dcfinitcncss, and expressed A t using (|10l) . 

Observe that new 77— dependent terms appear in the field equation for the scalar field compared to the scalar 
equation in 9(. They give a different z— dependent behavior of the scalar field, as we will discuss below. 

The asymptotic boundary condition (z — > 0) becomes 

*i(z) - Vi^ + V^ 2 ■ (18) 
At the horizon, we have the regularity condition 



2 - B{\ + i]B 2 ) 
(3-^)(l-^ 2 ) 



We will investigate the z— dependent profile of the scalar field perturbation, solving equation Q17p focusing 
on the case of scaling dimension A = 1 (setting ip% = 0). For a numerical solution, we impose the boundary 
condition = 1, and solve (jTTJ) with the use of a shooting method. Notice that the boundary condition at 

the horizon (|19j) depends on the coupling constant rj. The z— dependence of the scalar field for different values 
of the magnetic field is shown in Fig.[T] We see that as the magnetic field becomes stronger, the corresponding 
scalar perturbation in the bulk becomes bigger. This behavior is enhanced by the increasing strength of the 
interaction. 

In Fig. [2] we plot the phase diagram of the critical magnetic field and critical temperature. Above the lines, 
the holographic system will change into the normal phase. We use the critical temperature Tq for B = 0, rj = 
as a scale in the figure. The phase diagram reveals some very interesting features. 

Stronger magnetic fields corres pon d to lower critical temperatures. This behavior is expected, and for 
77 = it has been discussed in [lil llJ|. However in our case, as the strength of the higher-derivative coupling 
increases, smaller critical temperatures at strong enough magnetic field can be probed as can be seen on the 
right panel of Fig. [2l To see the behavior of the system at the quantum critical point (T = 0), we need to 
solve the full back-reacted system of Einstein-Maxwell-scalar field equations. Note that as the coupling 77 gets 



larger, for fixed critical temperature, the critical magnetic field B Cl above which the superconductor enters its 
normal phase, decreases. 

One of the reasons of introducing the generalized potential term ^ in condensed matter physics is to 
account for the FFLO states. These states appear in high-field superconductors when a strong magnetic 
field is applied, as Fulde and Ferrell 



17] showed independently. Their main 



la ] and Larkin and Ovchinnikov 
feature is exhibiting an order parameter which is not a constant, but it has a space variation. It was shown 
that if the coefficients in the functional ^ arc proportional to the magnetic field B, and the coefficient 7 of 
the gradient term in the Landau-Ginzburg functional is negative, then the maximum transition temperature 
corresponds to a space dependent state, and is larger than the transition temperature of the homogeneous 
state (for a detailed account see 5j). 

In our case, we recognize some similarities with the generalized Landau-Ginzburg theory. Looking at the 
right panel of Fig. [2] we see that in the low temperature limit and for a strong magnetic field our droplet 
solutions with negative coupling 77 have higher critical temperatures than the solutions with positive coupling. 

So far we concentrated on the strong magnetic field regime exhibited in Fig. [2] Let us now look at the left 
panel of Fig. [2] in the low magnetic field regime. We observe that, unlike strong magnetic fields, the critical 
temperature is larger if the higher-derivative coupling is stronger. This is consistent with the results obtained 
in section IIVI There is a transition region marked in the green box in the left part of Fig. [21 which we enlarge 
it in Fig. [3] We find that the cross point between lines of i] = —0.1 and i] = occurs at B ~ 0.8124. Below 
this value of the magnetic field, the critical temperature is smaller for negative coupling rj. Above this value, 
the critical temperature for ?; = —0.1 becomes higher than that of the minimal coupling. The cross point 
between lines of i] = —0.1 and i] = 0.1 appears at B « 0.8268. Above this transitional strength of the magnetic 
field, the critical temperature for r\ = —0.1 is higher than the positive coupling. The transitional strength 
of the magnetic field between 77 = 0.1 and the minimal coupling is at B « 0.8407, above which the critical 
temperature for minimal coupling becomes higher. 

Unlike in the Landau-Ginzburg theory 6j, we cannot compare the transition temperature of our inho- 
mogeneous solution to that of a homogeneous solution, because the latter does not exist in our system. 1 
Nevertheless, we can draw some conclusions, as we will discuss in the section IVl 



It is possible to compare homogeneous and inhomogeneous solutions in the presence of a magnetic field, if we introduce two 
gauge fields, as was discussed in [19]. The first gauge field couples to the scalar sourcing a charge condensate below a critical 
critical temperature, whereas the second gauge field incorporates a magnetic field that couples to the spin in the boundary 
theory. 
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FIG. 2: The critical temperature as a function of magnetic field. The low T behavior is enlarged in the right panel. 
The enlarged view in the green rectangle region is shown in Fig(3] 
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FIG. 3: The critical temperature as a function of magnetic field in the green rectangle region of Fig[2] 
IV. BELOW THE CRITICAL TEMPERATURE 



In this section we perform calculations below the critical temperature. We concentrate on the weak magnetic 
field limit. To simplify the calculations, we also focus on the probe limit by ignoring the back-reaction to the 
metric. Thus the metric takes the form of a planar Schwarzschild AdS black hole, 

1 



ds z 



dz 2 

-h{z)dt 2 + —— + dx 2 + dy 2 
h{z) 



h{z) = 1 



(20) 



We will study the condensation of the scalar field and calculate the critical temperature both analytically 
and numerically in the weak magnetic field limit A; — > (i = 1,2,3). This results in a solution which is 
approximately spatially independent. Thus, we consider a spherically symmetric ansatz 



* = $(z) , A„ = A t {z)dt . 



(21) 



Note that the potential term in the weak magnetic field limit (B —> 0) only contains gradients in the radial 
direction. Therefore, it merely induces interactions between the electric field and the scalar field. 



The Maxwell equation ([7]) and the scalar field equation (J5J in the black hole background (|2T)|) reduce to 



2r 1 z i ^A 2 {z^h' - 2h($ + zW))/rl 



2T 1 z 3 h^'(zh'^' + 2h(^' + zV")) 



-2rjhz 4 ^ 2 A 2 jr\ + z 2 h 2 (l + 2r]z 2 h^' 2 ) -2r]z^ 2 A 2 jr\ + z 2 h(l + 2r)z 2 h$>' 2 ) 



A' 



2V 2 (l + r 1 z i A' 2 /r 2 + ) 



-2r}z 4 ^ 2 A 2 /r 2 + + z 2 h(l + 2r 1 z 2 h^ /2 ) 



ti 2(1 + V z 4 A' 2 /rl + ^A'tA'l/rX) 
~h z(l- V z 4 A' 2 /r 2 + ) 



h 2 z 2 h(l - 7]z i Af/r\ 



A t = 



* = , 



(22) 



(23) 



where the prime denotes derivative in respect to z. Note that for vanishing coupling 77, the above field equations 



coincide with the corresponding equations in 1| 
infl 



As in [l| , the crucial feature is the last term in (|2"5|) which is a direct coupling of the scalar potential to the 
scalar field. This term comes from the effective mass term which reads 2 



z 2 A 2 {l-r 1 z i A' 2 /r 2 + ) 



r\h 



(24) 



and besides the scalar potential it also depends on the coupling constant 77. This is the gravitational analog 
7j of an Abelian U(l) symmetry breaking outside the horizon of the black hole background (|2"0]) compared to 
the symmetry breaking in the generalized Landau-Ginzburg theory of ([2]). 

To see the effect of the presence of the interaction term, in Fig. U]we plot the effective mass for m 2 = — 2 
and scaling dimension A = 1. The left panel shows that large negative values of the coupling constant r\ make 
m eff develop a deeper well. However, as we approach the horizon (right panel) larger values of the coupling 77 
give smaller values of the effective mass, which means that the condensation is easier to form. Therefore, the 
effect of the presence of the interaction term is that as the strength of the coupling 77 is increased, the system 
undergoes a phase transition in a higher critical temperature as can be seen in Table |U The same behavior 
can be observed in the case of scaling dimension A = 2 as can be seen in Fig. [5] and Table HU 

Next we will calculate the critical temperature analytically and show that it agrees with our numerical 
results. 

A. Analytic Results 



We will analytically determine the critical temperature for scaling dimension A = 1 with the use of the 
method developed in [8|. Since the scalar field vanishes at the critical temperature T c , the Maxwell equation 



2 m e f f i s coefficient of \l/ 2 in the action (SJ . 




FIG. 4: The value of the effective mass of the scalar field (|24[) as a function of z for scaling dimension A = 1 and 
various values of the coupling r\. 




FIG. 5: The value of the effective mass of the scalar field (|24[1 as a function of z for scaling dimension A = 2 and 
various values of the coupling 77. 



(f2"2"| reduces to A" — 0, which is solved by 



A t = Ar+(1 - z) , 



(25) 



where Arl is the charge density. On the other hand, near the boundary, we can introduce a function F(z) by 
defining 3 



* = zF(z) 

where F(z) has been normalized as F(Q) = 1 and F'(0) = . 

Thus, substituting (|2"5"|) and (|26p into the scalar field equation (f2"5)l we obtain 



(26) 



F"(z) 



5z 3 -2 2(1 + 77AV) 



z(l-z 3 ) z(l-r 1 \ 2 z i ) 
A 2 



F\z)- 



3z 



2(1 + t 1 X 2 z 4 ) 



1-z 3 ^(l-^z 4 ) Z 2 (1-Z 3 )(1-7?Z 4 A 2 ) 



(1 + z + z 3 ) 2 



F = 



F(z) 
(27) 



For the other scaling dimension, A = 2, the definition becomes if = ^2 2 ^ z 2 F(z). 



The above equation can be converted into 
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[T(z)F>(z)}'-Q(z)F(z) = 



(28) 



with 



T{z) = (1 - z A ){\ - r/X z ) , Q(z) = z - A 



»(l-z)(l-i,AV) 



?7AV(5,z 3 -2) 



1 + 2 + Z 2 

With the use of Sturm-Liouville theory, the eigenvalue A can be found as the minimum solution of 



[ dz [T{z)F'{z) 2 + Q(z)F{zf] =0 
Jo 



To proceed, we assume the trial function 



(29) 



(30) 



F(z) ee 1 - az 2 , 



(31) 



which satisfies the condition F(0) = 1 as well as F'(0) = 0. After integrating over z in (f30|) . we find that A 2 
satisfies the second-order equation 



02 A 4 + a\\ 2 + ao = 



(32) 



where 



a-2 = V < \ a 



ai = rja 



409 

2V3tt + 3 ( -— + 21n3 



2a 



2V3ir - 3 ( j -f . 



13- 12 In 3 



-rja 



110501 
2310 
3263 



105 



f &V3tt + 18 In 3 
12V3tt - 36 In 3 



- a 2 (12 In 3 - 13) + a(4V37r + 12 In 3 - 36) 
3r?fy-121n3j +2V37r-61n3 



T0a 2 + 6a-6 . 



(33) 



For a fixed value of the coupling 77, we obtain the minimum value of A 2 by varying a. Then the critical 
temperature is found to be 



(34) 



B. Numerical Results 



The field equations (|22|) and (|23|) were solved numerically with m 2 = — 2. In Table U we list the critical 
temperature calculated numerically for various values of the coupling constant 77, and compare it with the 
value of the critical temperature calculated using the analytic result ([34|) . Note that for r\ = we recover 



12 




T] =-0.3 
T] =-0.01 
r,=0 
i] =0.01 

T) =0.03 




FIG. 6: The order parameters (d) (left panel) and (O2) (right panel) as functions of temperature for various values 
of the coupling 77. 



n 


-0.3 


-0.2 


-0.1 





0.01 


0.1 


0.2 


0.3 


T c /y/~p (numerical) 


0.2224 


0.2234 


0.2244 


0.2255 


0.2257 


0.2267 


0.2280 


0.2294 


Tel \fp (analytic) 
(a) 


0.2218 
(0.2868) 


0.2228 
(0.2720) 


0.2239 
(0.2561) 


0.2250 
(0.2389) 


0.2251 
(0.2371) 


0.2261 
(0.2204) 


0.2273 
(0.2006) 


0.2287 
(0.1793) 



TABLE I: The critical temperature for scaling dimension A = 1 and various values of the coupling 77. 

the results of |l|. We also list the corresponding values of a. From this table, we can see that the numerical 
results and the analytic results are in good agreement. In Table fill we also summarize the critical temperature 
for the case of scaling dimension A = 2 calculated analytically. We see that the critical temperature increases 
as the coupling 77 increases. 

Next we examine the effect of the presence of the higher-derivative coupling term on the condensation of the 
scalar field. In Fig.[6]we show the vacuum expectation values of the two operators 0\ and O2 versus the critical 
temperature. We see that with the increase of the strength 77 of the interaction, the gap becomes smaller, 
which means that the scalar operator can condense easier when the coupling 77 is stronger. This agrees with 
the property of the critical temperature we discussed above. For higher 77, our numerical calculation becomes 
harder at lower temperature, however this does not prevent us from observing the qualitative influence of the 
interaction on the condensation. For scaling dimension A = 1, this behavior is more pronounced (see left 
panel of Fig. [6j . For scaling dimension A = 2, we enlarge a local region inserted into the right panel of Fig. |6] 
to bring out this behavior. From the critical temperature behavior obtained with both numerical and analytic 
methods, together with the observed gap behavior, we can conclude that in the weak magnetic field limit, the 
greater the strength of the interaction the easier it is for the condensation to form. 





-0.2 


-0.1 





0.01 


0.1 


0.2 


0.3 


Tc/yfp (analytic) 


0.1092 


0.1157 


0.1170 


0.1170 


0.1191 


0.1384 


0.1506 



TABLE II: The critical temperature for scaling dimension A = 2 and various values of the coupling 77. 
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V. DISCUSSION 

Motivated by the generalized Landau-Ginzbourg theory, we set up a gravitational dual of a holographic 
superconductor which included a higher-derivative coupling of strength r\. Our system had standard compo- 
sition, consisting of a f (1) gauge field, and a complex scalar field coupled to gravity. The novel feature was 
a higher-derivative coupling between the U(l) gauge field and the scalar. We solved the field equations and 
compared with the conventional case of a holographic superconductor (rj = 0). 

In the limit in which the magnetic field is weak, we found a spatially independent (homogeneous) solution. 
In the case in which the coupling r\ of the higher-derivative coupling vanishes, we recovered the results of 
Q. As the strength of the coupling increases, the gravitational mechanism of breaking an Abelian U(l) 
symmetry outside the horizon of a black hole [7j becomes more effective, as can be seen in Fig. 2] and Fig. [5] 
Calculating the critical temperature, both analytically and numerically, we found that as the strength of the 
interaction increased, there was an enhancement of the critical temperature (Tables U and [TTJ) at which the 
scalar condensate formed. 

In the presence of a magnetic field, the normal state corresponds to the background of a dyonic black hole. 
The field equations possess spatially dependent (inhomogeneous) solutions. We found that the presence of 
the higher-derivative coupling did not affect the x-dependent profile of the scalar field solution. The evolution 
of the x-dependent profile of the scalar field is given by the solution of equation ([16]) which depends on the 
strength of the magnetic field. As it is discussed in 9|, when the magnetic field is zero, the x-dependence of the 
condensate disappears and we recover the homogeneous solution. However, as the magnetic field is increased, 
the cc-depcndcnt profile of the scalar field acquires a "thickness" making the condensate more inhomogeneous. 

The presence of the higher-derivative coupling term influences the radial dependence of the scalar field 
contributing to the critical temperature. In Fig. [2 the (B, T) phase diagram for various values of the coupling 
r\ is shown. We can see some interesting features of our droplet solutions in the low temperature, strong 
magnetic field region of the phase diagram (right panel of Fig. [2]). For 77 > and a fixed value of the magnetic 
field we can probe lower temperatures. For 77 < the condensate forms at a higher critical temperature 
and at a stronger magnetic field making the corresponding droplet solution more inhomogeneous. As we 
discussed, this behavior is reminiscent of the formation of FFLO states in the Landau-Ginzburg theory, where 
the inhomogeneous states at low temperature are energetically more favorable. 

To have a better understanding of our droplet solutions in the T — > limit and their possible connection 
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with the inhomogeneous FFLO states, we have to improve our numerical techniques in order to probe larger 
(positive and negative) values of the coupling ?/. We could also look for solutions with a strong magnetic field 
penetrating the whole a?y-plane and study the behavior of the system by solving the full non-linear system of 
the coupled Einstein-Maxwell-scalar equations. 
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Abstract 

We discuss the gravitational dual of a holographic superconductor consisting of a (7(1) gauge field, 
a complex scalar field coupled to a charged AdS black hole and a higher-derivative coupling between 
the U(l) gauge field and the scalar with coupling constant 77. In the presence of a magnetic field, the 
system possesses localized spatially dependent droplet solutions which, in the low temperature limit, 
have smaller critical temperature for 77 > than the droplet solutions without the interaction term 
(rj — 0). In the weak magnetic field limit, the opposite behavior is observed: the critical temperature 
increases as we increase n. We also calculate the energy gap in the probe limit and find that it is 
larger for 77 < than the energy gap in the conventional case (n — 0). 
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I. INTRODUCTION 



AO 



The plicnomonological holographic superconductor [1|, |2| has many similarities with the effective Landau- 
Ginzburg theory. The gravity sector of the holographic superconductor is an Abelian-Higgs model with a 
dynamical black hole metric and the scalar field condenses below a certain critical temperature. The Landau- 
Ginzburg theory is the effective field theory description, in one less dimension, of a superconductor near the 
superconducting phase transition. The dynamical degree of freedom is the complex order parameter, ^ = (0), 
which is coupled to a background electromagnetic field. This is described by the standard Ginzburg-Landau 
functional 

F = a\^ + ^\\ (1) 

and the coefficient a vanishes at the transition temperature T = T c . At T < T c , the coefficient a is negative 
and the minimum of F occurs for a uniform superconducting state at |^| 2 = — ^. 

However, the Landau-Ginzburg theory is not a microscopic theory and for the description of the dynamics 
of a Cooper pair the BCS theory is required. By the same token, one hopes that the phenomenological 
description of holographic superconductors will find its microscopic realization in string theory 3|, m- I n 
spite of this limitation, if one insists on working with the Landau-Ginzburg theory and wishes to describe 
superconductivity away from the critical temperature and in the presence of a magnetic field.the functional 
([T]) has to be supplemented with higher order terms involving spatial derivatives of the form 5| 

Fq = a(B, T)\V\ 2 + 7 (g, T)|^| 2 + ^ T) |^ 2 fr| 2 + ^ T) |^| 4 . (2) 

Note that the various coefficients also depend on the magnetic field. The critical temperature resulting from 
the generalized functional (|2|) depends on the value of the coefficient 7 and for 7 < interesting new physics 
emerges 6 1 . 

Despite the similarities between the holographic superconductors and the Landau-Ginzburg theory, there 
are important differences discussed in [2(. However, to keep the similarities of the holographic superconductor 
with the generalized Landau-Ginzburg theory, we can add a potential term which depends on the scalar field 
'J, as well as its spatial derivatives along the boundary, 

y(*,^ 2 *,V 2 *,...) . (3) 



While spatial derivatives are used in @, they are problematic in our case, because they signal the breakdown 
of Lorentz invariance. To obtain an expression similar to @ in a covariant way, we note that the vector 
F^ v d u ^ consists of the derivatives of ^> along the x and y directions, if the electromagnetic field strength F M „ 
describes a magnetic field in the radial direction (perpendicular to the boundary) . This leads us to consider 



(4) 



where F = dA is the field strength of a U(l) gauge field A^, ^ is a charged complex scalar field of charge q 
and mass m and = — iqA^. 

In this work we study the effects of the presence of the higher-derivative terms on the holographic super- 
conductor. To avoid technical difficulties, we set all couplings, except r\ to zero (rf = ■ ■ ■ = 0). In section 
HTl we present the details of the model. In section IIII1 we consider droplet solutions and calculate the critical 
temperature. In section IIV1 we go below the critical temperature and calculate the energy gap in the probe 
limit. Finally in section |V1 we discuss our conclusions. 



II. THE MODEL 



We consider the action 



S = d 4 c 



R + 6/L 2 1 



-F^F^ - \D^\ 2 - VpS!) 



16ttG 

where the potential term is given by ^ with rj' = ■ ■ ■ = 0. The field equations are: 
• the Einstein equations 



R 



1 3 



(5) 



(6) 



• the Maxwell equations 



V^F^ + ^=<9 M [y/^g{D K ^f)(D x ^y (g KV F» x - g K »F uX + g^F^ - g^ x F UK )] 
= iq - + iqr]g iip F pv [F^*(D K ^) - F^ A *(L> A *)*] 



(7) 



• and the scalar field equation 
1 



-9 



d ll [y/=$g> iv {d v V - iqA v *)] + iqgT ' A v (d^ - iqA^) + m 2 4- 



V 



■A [y/=g~g KX F KV F x » {d„* - iqA v ^)] + iqr,g KX F KV F x v A„ (d v V - iqA,^) = . (8) 



In this work, we will set L = 1, 8irG = 1, g = 1. 

We note that the presence of the coupling constant r\ adds new terms in the field equations which makes 
the system of the differential equations highly non-trivial. However, in spite of the complexity of the system, 
we were able to perform numerical calculations and also calculate the critical temperature analytically in the 
weak field limit. 

III. THE CRITICAL TEMPERATURE 



The original motivation of introducing a potential term of the form ^ in the Landau-Ginzgurg theory was 
to study the behavior of a superconductor in the presence of a strong magnetic field and the generation of 
anisotropic states in the low temperature limit. To investigate further similarities of a holographic supercon- 
ductor to the generalized Landau-Ginzburg theory, in this section we switch on a magnetic field and study 
the effect of spatially dependent solutions on the critical temperature. 



Spatially dependent solutions in the presence of a magnetic field were studied in [9h11| . Solving the Maxwell 



scalar equations in a dyonic black hole background localized droplet and vortex solutions were found and their 



(B,T) phase diagram was studied. In 



with the addition of an external magnetic field, it was argued 



that the holographic superconductor is of type II as the magnetic field is lowered. An analytic study of a 



holographic superconductor was carried out in 



12j . where the upper critical magnetic field was calculated 



analytically. Vortex lattice solutions in the presence of a magnetic field were presented in 13 1. 

To study the effect of the interaction term (with r\ ^ 0), we will consider the backreaction of the Maxwell 
field to gravity. The Einstein-Maxwell field equations ((6]) and (0 at the critical temperature are solved by a 



dyonic black hole (i.e., one with both an electric and a magnetic charge). The metric reads 



14] 



ds 2 



dr 2 



+ r 2 [-h{r)dt 2 + dx 2 + dy 2 



h(r] 



(X 2 + B 2 )r% (I 
4r 3 



The gauge field is given by 



Aft = A t {r)dt + A y (x)dy = Ar+ ( 1 - ^ ) dt + Br\xdy , 



(9) 



(10) 

where according to the AdS / CFT dictionary, the charge density and the chemical potential on the boundary 
are given by p = Ar^_ and /x = Ar + , respectively. Br 2 ^ is the magnitude of the external magnetic field (in the 
radial direction, r, perpendicular to the boundary xy-plane). At the horizon, r — r+, we have At = 0. The 
Hawking temperature of the black hole is 



3r + / A 2 +B 2 
4tt I 12 



(11) 



In the above background with a magnetic field B, we will investigate the spatially dependent solutions of the 
field equations by considering the ansatz ^ = ^(r, x). Near the AdS boundary, the asymptotic behavior is 
^ = -^L- + with A+ = | ± i-v/9 + 4m 2 . Both and "f? correspond to normalizablc modes, therefore we 

r — T ' J, Z 

obtain two different systems (labeled by A±) in which one mode is the source and the other gives the vacuum 
expectation value of the dual operator of dimension Aj ('J'i ~ (Cj))- 

Right below the critical temperature, we may regard the scalar field as a perturbation in the gravitational 
background of © and (|10l) . Then the scalar field equation in the coordinate z 



1 + 

h 













-!)■ 







,4 /i/2 iJ2 



becomes 



(1 V z*A?/t2)Al gV ^ v 



ft 



(12) 



where the prime denotes differentiation with respect to z. Without loss of generality, we can take to be real 
and separate 'J into functions of a single variable, 



Combining equations (fT2]) and (p~3|) . we obtain 
ft 



(13) 



ft' 2 
"ft ~ "z 



ft*; 



(i-')^; 2 K)^ 



?*ift 



^(l+^V) 

*2 L J 



For consistency of this equation, we need 



#2 - £ 2 z 2 *2 = -fc 2 ^ 



(14) 



(15) 



where k 2 is an arbitrary parameter. This equation is of the form of Schrodinger's equation for a simple 
harmonic oscillator. The eigenvalues are 



fc 2 = (2n+ l)B , n = 0,1,2,. 



(16) 



and the corresponding eigenfunctions are given in terms of Hermite polynomials, H n (y2Bx). Note that the 
x-dependent part of the scalar field obeys the same equation as the one found in 9|. As explained in {J, we 
ought to choose the ground state (n = 0) to find the critical temperature. 

Substituting the x-dependent equation (JT5J) into (|T4)) . we obtain the field equation for ty\{z) which reads 

'ft' 2 



[1 - V X 2 z 4 ] *i' 
2 

j2ft 



2 „4 



r/X z 



ft' 2 

ft z y ' \ h z 
X 2 (l-i 1 X 2 z 4 )(l-z) 2 _ 5(1 + ?? ff 2 z 4 ) 
ft 2 ft 



*i = , 



(17) 



FIG. 1: The ^-dependence of the scalar field for various values of the magnetic field B and coupling constant rj. 

where we have chosen m 2 = — 2 for dcfinitcncss, and expressed A t using (|10l) . 

Observe that new 77— dependent terms appear in the field equation for the scalar field compared to the scalar 
equation in 9(. They give a different z— dependent behavior of the scalar field, as we will discuss below. 

The asymptotic boundary condition (z — > 0) becomes 

*i(z) - Vi^ + V^ 2 ■ (18) 
At the horizon, we have the regularity condition 



2 - B{\ + i]B 2 ) 
(3-^)(l-^ 2 ) 



We will investigate the z— dependent profile of the scalar field perturbation, solving equation Q17p focusing 
on the case of scaling dimension A = 1 (setting ip% = 0). For a numerical solution, we impose the boundary 
condition = 1, and solve (jTTJ) with the use of a shooting method. Notice that the boundary condition at 

the horizon (|19j) depends on the coupling constant rj. The z— dependence of the scalar field for different values 
of the magnetic field is shown in Fig.[T] We see that as the magnetic field becomes stronger, the corresponding 
scalar perturbation in the bulk becomes bigger. This behavior is enhanced by the increasing strength of the 
interaction. 

In Fig. [2] we plot the phase diagram of the critical magnetic field and critical temperature. Above the lines, 
the holographic system will change into the normal phase. We use the critical temperature Tq for B = 0, rj = 
as a scale in the figure. The phase diagram reveals some very interesting features. 

Stronger magnetic fields corres pon d to lower critical temperatures. This behavior is expected, and for 
77 = it has been discussed in [lil llJ|. However in our case, as the strength of the higher-derivative coupling 
increases, smaller critical temperatures at strong enough magnetic field can be probed, as can be seen on the 
right panel of Fig. [2 To see the behavior of the system at the quantum critical point (T = 0), we need to 
solve the full back-reacted system of Einstein-Maxwell-scalar field equations. Note that as the coupling 77 gets 



larger, for fixed critical temperature, the critical magnetic field B Cl above which the superconductor enters its 
normal phase, decreases. 

One of the reasons of introducing the generalized potential term ^ in condensed matter physics is to 
account for the FFLO states. These states appear in high-field superconductors when a strong magnetic 
field is applied, as Fulde and Ferrell 



17] showed independently. Their main 



la ] and Larkin and Ovchinnikov 
feature is exhibiting an order parameter which is not a constant, but it has a space variation. It was shown 
that if the coefficients in the functional ^ arc proportional to the magnetic field B, and the coefficient 7 of 
the gradient term in the Landau-Ginzburg functional is negative, then the maximum transition temperature 
corresponds to a space dependent state, and is larger than the transition temperature of the homogeneous 
state (for a detailed account see 5j). 

In our case, we recognize some similarities with the generalized Landau-Ginzburg theory. Looking at the 
right panel of Fig. [2] we see that in the low temperature limit and for a strong magnetic field our droplet 
solutions with negative coupling 77 have higher critical temperatures than the solutions with positive coupling. 

So far we concentrated on the strong magnetic field regime exhibited in Fig. [2] Let us now look at the left 
panel of Fig. [2] in the low magnetic field regime. We observe that, unlike strong magnetic fields, the critical 
temperature is larger if the higher-derivative coupling is stronger. This is consistent with the results obtained 
in section IIVI There is a transition region marked in the green box in the left part of Fig. [21 which we enlarge 
it in Fig. [3j We find that the cross point between lines of 77 = —0.1 and 77 = occurs at B/Tq sa 21.74. Below 
this value of the magnetic field, the critical temperature is smaller for negative coupling 77. Above this value, 
the critical temperature for 77 = —0.1 becomes higher than that of the minimal coupling. The cross point 
between lines of 77 = —0.1 and 77 = 0.1 appears at B/Tq sa 21.91. Above this transitional strength of the 
magnetic field, the critical temperature for 77 = —0.1 is higher than the positive coupling. The transitional 
strength of the magnetic field between 77 = 0.1 and the minimal coupling is at B/Tq rj 22.15, above which the 
critical temperature for minimal coupling becomes higher. 

Unlike in the Landau-Ginzburg theory |6[, we cannot compare the transition temperature of our inho- 
mogeneous solution to that of a homogeneous solution, because the latter does not exist in our system. 1 
Nevertheless, we can draw some conclusions, as we will discuss in the section IVl 



It is possible to compare homogeneous and inhomogeneous solutions in the presence of a magnetic field, if we introduce two 
gauge fields, as was discussed in [19]. The first gauge field couples to the scalar sourcing a charge condensate below a critical 
critical temperature, whereas the second gauge field incorporates a magnetic field that couples to the spin in the boundary 
theory. 
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FIG. 2: The critical temperature as a function of magnetic field. The low T behavior is enlarged in the right panel. 
The enlarged view in the green rectangle region is shown in Fig(3] 
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FIG. 3: The critical temperature as a function of magnetic field in the green rectangle region of Fig[2] 
IV. BELOW THE CRITICAL TEMPERATURE 



In this section we perform calculations below the critical temperature. We concentrate on the weak magnetic 
field limit. To simplify the calculations, we also focus on the probe limit by ignoring the back-reaction to the 
metric. Thus the metric takes the form of a planar Schwarzschild AdS black hole, 

1 



ds z 



dz 2 

-h{z)dt 2 + —— + dx 2 + dy 2 
h{z) 



h{z) = 1 



(20) 



We will study the condensation of the scalar field and calculate the critical temperature both analytically 
and numerically in the weak magnetic field limit A4 — > (i = 1,2,3). This results in a solution which is 
approximately spatially independent. Thus, we consider a spherically symmetric ansatz 



* = $(z) , A„ = A t {z)dt . 



(21) 



Note that the potential term in the weak magnetic field limit (B —> 0) only contains gradients in the radial 
direction. Therefore, it merely induces interactions between the electric field and the scalar field. 



The Maxwell equation ([7]) and the scalar field equation (J5J in the black hole background (|2T)|) reduce to 



2r 1 z i ^A 2 {z^h' - 2h($ + zW))/rl 



2r]z 3 h$>'(zh , ty' + 2h{^' + zW))/r\ 



-2rjhz 4 ^ 2 A 2 jr\ + z 2 h 2 (l + 2r]z 2 h^' 2 ) -2r]z^ 2 A 2 /r 2 + + z 2 h(l + 2r)z 2 h$>' 2 ) 



A' 



2V 2 (l + r 1 z i A' 2 /r 2 + ) 



-2r}z 4 ^ 2 A 2 /r 2 + + z 2 h(l + 2r 1 z 2 h^ /2 ) 



ti 2(1 + V z 4 A' 2 /rl + ^A'tA'l/rX) 
~h z(l- V z 4 A' 2 /r 2 + ) 



A} 



h 2 z 2 h(l - 7]z i Af/r\ 



A t = 



* = , 



(22) 



(23) 



where the prime denotes derivative in respect to z. Note that for vanishing coupling 77, the above field equations 



coincide with the corresponding equations in 1| 
infl 



As in [l| , the crucial feature is the last term in (|2"5|) which is a direct coupling of the scalar potential to the 
scalar field. This term comes from the effective mass term which reads 2 



z 2 A 2 {l-r 1 z i A' 2 /r 2 + ) 



r\h 



(24) 



and besides the scalar potential it also depends on the coupling constant 77. This is the gravitational analog 
7j of an Abelian U(l) symmetry breaking outside the horizon of the black hole background (|2"0]) compared to 
the symmetry breaking in the generalized Landau-Ginzburg theory of ([2]). 

To see the effect of the presence of the interaction term, in Fig. U]we plot the effective mass for m 2 = — 2 
and scaling dimension A = 1. The left panel shows that large negative values of the coupling constant r\ make 
m eff develop a deeper well. However, as we approach the horizon (right panel) larger values of the coupling 77 
give smaller values of the effective mass, which means that the condensation is easier to form. Therefore, the 
effect of the presence of the interaction term is that as the strength of the coupling 77 is increased, the system 
undergoes a phase transition in a higher critical temperature as can be seen in Table |U The same behavior 
can be observed in the case of scaling dimension A = 2 as can be seen in Fig. [5] and Table HU 

Next we will calculate the critical temperature analytically and show that it agrees with our numerical 
results. 

A. Analytic Results 



We will analytically determine the critical temperature for scaling dimension A = 1 with the use of the 
method developed in [8|. Since the scalar field vanishes at the critical temperature T c , the Maxwell equation 



2 m e f f i s coefficient of \l/ 2 in the action (SJ . 




FIG. 4: The value of the effective mass of the scalar field (|24[) as a function of z for scaling dimension A = 1 and 
various values of the coupling r\. 




FIG. 5: The value of the effective mass of the scalar field (|24[1 as a function of z for scaling dimension A = 2 and 
various values of the coupling 77. 



(f2"2"| reduces to A" — 0, which is solved by 



A t = Ar+(1 - z) , 



(25) 



where Arl is the charge density. On the other hand, near the boundary, we can introduce a function F(z) by 
defining 3 



* = zF(z) 

where F(z) has been normalized as F(Q) = 1 and F'(0) = . 

Thus, substituting (|2"5"|) and (|26p into the scalar field equation (f2"5)l we obtain 



(26) 



F"(z) 



5z 3 -2 2(1 + 77AV) 



z(l-z 3 ) z(l-r 1 \ 2 z i ) 
A 2 



F\z)- 



3z 



2(1 + t 1 X 2 z 4 ) 



1-z 3 ^(l-^z 4 ) Z 2 (1-Z 3 )(1-7?Z 4 A 2 ) 



(1 + z + z 3 ) 2 



F = 



F(z) 
(27) 



For the other scaling dimension, A = 2, the definition becomes if = ^2 2 ^ z 2 F(z). 



The above equation can be converted into 



11 



[T(z)F>(z)}'-Q(z)F(z) = 



(28) 



with 



T{z) = (1 - z A ){\ - r/X z ) , Q(z) = z - A 



»(l-z)(l-i,AV) 



?7AV(5,z 3 -2) 



1 + 2 + Z 2 

With the use of Sturm-Liouville theory, the eigenvalue A can be found as the minimum solution of 



[ dz [T{z)F'{z) 2 + Q(z)F{zf] =0 
Jo 



To proceed, we assume the trial function 



(29) 



(30) 



F(z) ee 1 - az 2 , 



(31) 



which satisfies the condition F(0) = 1 as well as F'(0) = 0. After integrating over z in (f30|) . we find that A 2 
satisfies the second-order equation 



02 A 4 + a\\ 2 + ao = 



(32) 



where 



a-2 = V < \ a 



ai = rja 



409 

2V3tt + 3 ( -— + 21n3 



2a 



2V3ir - 3 ( j -f . 



13- 12 In 3 



-rja 



110501 
2310 
3263 



105 



f &V3tt + 18 In 3 
12V3tt - 36 In 3 



- a 2 (12 In 3 - 13) + a(4V37r + 12 In 3 - 36) 
3r?fy-121n3j +2V37r-61n3 



T0a 2 + 6a-6 . 



(33) 



For a fixed value of the coupling 77, we obtain the minimum value of A 2 by varying a. Then the critical 
temperature is found to be 



(34) 



B. Numerical Results 



The field equations (|22|) and (|23|) were solved numerically with m 2 = — 2. In Table U we list the critical 
temperature calculated numerically for various values of the coupling constant 77, and compare it with the 
value of the critical temperature calculated using the analytic result ([34|) . Note that for r\ = we recover 
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r,=0 
i] =0.01 

T) =0.03 




FIG. 6: The order parameters (d) (left panel) and (O2) (right panel) as functions of temperature for various values 
of the coupling 77. 



n 


-0.3 


-0.2 


-0.1 





0.01 


0.1 


0.2 


0.3 


T c /y/~p (numerical) 


0.2224 


0.2234 


0.2244 


0.2255 


0.2257 


0.2267 


0.2280 


0.2294 


Tel \fp (analytic) 
(a) 


0.2218 
(0.2868) 


0.2228 
(0.2720) 


0.2239 
(0.2561) 


0.2250 
(0.2389) 


0.2251 
(0.2371) 


0.2261 
(0.2204) 


0.2273 
(0.2006) 


0.2287 
(0.1793) 



TABLE I: The critical temperature for scaling dimension A = 1 and various values of the coupling 77. 

the results of |l|. We also list the corresponding values of a. From this table, we can see that the numerical 
results and the analytic results are in good agreement. In Table fill we also summarize the critical temperature 
for the case of scaling dimension A = 2 calculated analytically. We see that the critical temperature increases 
as the coupling 77 increases. 

Next we examine the effect of the presence of the higher-derivative coupling term on the condensation of the 
scalar field. In Fig.[6]we show the vacuum expectation values of the two operators 0\ and O2 versus the critical 
temperature. We see that with the increase of the strength 77 of the interaction, the gap becomes smaller, 
which means that the scalar operator can condense easier when the coupling 77 is stronger. This agrees with 
the property of the critical temperature we discussed above. For higher 77, our numerical calculation becomes 
harder at lower temperature, however this does not prevent us from observing the qualitative influence of the 
interaction on the condensation. For scaling dimension A = 1, this behavior is more pronounced (see left 
panel of Fig. [6j . For scaling dimension A = 2, we enlarge a local region inserted into the right panel of Fig. |6] 
to bring out this behavior. From the critical temperature behavior obtained with both numerical and analytic 
methods, together with the observed gap behavior, we can conclude that in the weak magnetic field limit, the 
greater the strength of the interaction the easier it is for the condensation to form. 





-0.2 


-0.1 





0.01 


0.1 


0.2 


0.3 


Tc/yfp (analytic) 


0.1092 


0.1157 


0.1170 


0.1170 


0.1191 


0.1384 


0.1506 



TABLE II: The critical temperature for scaling dimension A = 2 and various values of the coupling 77. 
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V. DISCUSSION 

Motivated by the generalized Landau-Ginzbourg theory, we set up a gravitational dual of a holographic 
superconductor which included a higher-derivative coupling of strength r\. Our system had standard compo- 
sition, consisting of a f (1) gauge field, and a complex scalar field coupled to gravity. The novel feature was 
a higher-derivative coupling between the U(l) gauge field and the scalar. We solved the field equations and 
compared with the conventional case of a holographic superconductor (rj = 0). 

In the limit in which the magnetic field is weak, we found a spatially independent (homogeneous) solution. 
In the case in which the coupling r\ of the higher-derivative coupling vanishes, we recovered the results of 
Q. As the strength of the coupling increases, the gravitational mechanism of breaking an Abelian U(l) 
symmetry outside the horizon of a black hole [7j becomes more effective, as can be seen in Fig. 2] and Fig. [5] 
Calculating the critical temperature, both analytically and numerically, we found that as the strength of the 
interaction increased, there was an enhancement of the critical temperature (Tables U and [TTJ) at which the 
scalar condensate formed. 

In the presence of a magnetic field, the normal state corresponds to the background of a dyonic black hole. 
The field equations possess spatially dependent (inhomogeneous) solutions. We found that the presence of 
the higher-derivative coupling did not affect the x-dependent profile of the scalar field solution. The evolution 
of the x-dependent profile of the scalar field is given by the solution of equation ([16]) which depends on the 
strength of the magnetic field. As it is discussed in 9|, when the magnetic field is zero, the x-dependence of the 
condensate disappears and we recover the homogeneous solution. However, as the magnetic field is increased, 
the cc-depcndcnt profile of the scalar field acquires a "thickness" making the condensate more inhomogeneous. 

The presence of the higher-derivative coupling term influences the radial dependence of the scalar field 
contributing to the critical temperature. In Fig. [2 the (B, T) phase diagram for various values of the coupling 
r\ is shown. We can see some interesting features of our droplet solutions in the low temperature, strong 
magnetic field region of the phase diagram (right panel of Fig. [2]). For 77 > and a fixed value of the magnetic 
field we can probe lower temperatures. For 77 < the condensate forms at a higher critical temperature 
and at a stronger magnetic field making the corresponding droplet solution more inhomogeneous. As we 
discussed, this behavior is reminiscent of the formation of FFLO states in the Landau-Ginzburg theory, where 
the inhomogeneous states at low temperature are energetically more favorable. 

To have a better understanding of our droplet solutions in the T — > limit and their possible connection 
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with the inhomogeneous FFLO states, we have to improve our numerical techniques in order to probe larger 
(positive and negative) values of the coupling ?/. We could also look for solutions with a strong magnetic field 
penetrating the whole a?y-plane and study the behavior of the system by solving the full non-linear system of 
the coupled Einstein-Maxwell-scalar equations. 
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